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Computation  of  the 

Residual  Radiation  Transmission  Factor  for 
the  Marine  Corps  LVTP7 


by 

Robert  G.  Cawley 


We  develope  an  analytical  technique  for  computing  the  gamma 
radiation  reduction  factor  for  a  slab  of  arbitrary  orientation  in  a 
field  of  uniformly  distributed  radioactive  fallout  (Co60 ) .  We 
apply  this  to  a  model  of  the  army  M113  APC,  assuming  the  entire 
reduction  is  due  to  material  in  the  vehicle  walls,  ceiling  and  floor. 
The  r^cnlting  protection,  P  =  1.5,  differs  sharply  with  the  measured 
value,  P  =  4.0,  which  we  assume  is  due  to  the  presence  of  other 
material  than  that  in  the  vehicle  surfaces.  As  a  second  model  of 
the  M113  we  assume  its  shape  to  be  that  of  a  parallelepiped  having 
a  square  vertical  cross  section  of  side  L,  uniform  wall  thickness 
and  height  equal  to  the  actual  vehicle  height.  We  choose  L 

and  t  so  that  the  surface  area  and  total  mass  ->f  the  model  vehicle 

ef  r 

agree  with  those  of  the  actual  vehicle.  A  10X  albedo  correction  to 
the  protection  then  results  in  agreement  with  experiment: 

^theory  =  ^  ^ *  We  aPPly  the  same  approach  to  the  LVTP7,  getting 
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^theory  =  4.1.  The  radiation  reduction  factors,  R  =  P,  for  these 

two  cases  are  0.26  and  0.24  respectively;  and  the  residual  radiation 
transmission  factors,  defined  from  R  to  include  the  effects  of  a 
ground  roughness  parameter  of  3  ft,  are  both  0.3. 
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Computation  of  the  Residual  Radiation  Transmission  Factor  for  the 
Marine  Corps  LVTP7. 

This  report  gives  a  new  method  for  calculating  the  y-radiation 
transmission  through  a  plane  slab  of  arbitrary  orientation,  above  a 
field  of  uniformly  distributed  radioactive  fallout.  The  method  is 
applied  to  the  Marine  Corps  LVTP7  to  estimate  the  residual  radiation 
transmission  factor.  The  work  was  sponsored  by  the  U.  S.  Marine 
Corps  Development  and  Education  Command. 

The  author  wishes  to  thank  M.  A.  Schmoke  for  a  discussion,  and  for 
securing  certain  data  pertaining  to  the  army  M113  APC,  and  C.  M. 
Huddleston  for  discussions  and  for  assistance  in  obtaining  reference 
material,  and  for  reading  the  report. 
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I.  INTRODUCTION 

The  present  work  reports  the  results  of  an  effort  to  compute  the 
protection  factor  and  the  residual  radiation  transmission  factor 
(RRTFrdefined  in  Section  2)  for  the  Marine  Corps  LVTP7  troop  carrier 
in  a  field  of  radioactive  fallout.  The  P7  is  similar  to  another 
vehicle,  the  army  M113  APC  {armored  personnel  carrier),  for  which 
experimental  data  are  available,1  and  we  have  used  those  results  in 
our  method.  Briefly,  we  have  developed  an  approximate  analytical 
technique,  which  enables  us  to  calculate  the  protection  factor  for 
an  enclosure  whose  walls  are  composed  of  plane  slabs  having  arbitrary 
orientations. 

Applied  to  a  fairly  faithful  model  of  the  M113  geometry  and  using 
for  wall  thicknesses  the  actual  values,  the  calculated  value  of  the 
protection  factor  is  P  =  1.53,  much  lower  than  the  experimental 
value  of  4.0.  The  reason  for  this  discrepancy  is  that  the  mass  of 
the  M113  in  the  model  is  less  than  half  the  mass  of  the  actual 
vehicle  in  the  {combat  loaded)  experimental  conditions,  so  there  is 
quite  a  bit  more  "extra"  material  actually  present  to  enhance  the 
protection  over  that  afforded  by  the  model;  this  material  includes 
the  engine,  treads,  wheels,  axles,  fuel,  troops,  interior  panelling 
surfaces,  and  sundry  paraphernalia,  the  mass  distribution  of  which 
it  is  impractical  to  try  to  take  into  account.  In  the  first  place 
the  labor  involved  would  be  enormous,  and  the  results  withal 
probably  not  trustworthy.  Even  more  significantly,  there  is  a 
"theoretical  error"  inherent  in  the  notion  of  vehicle  protection 
factor,  a  circumstance  which  makes  too  precise  a  result  meaningless: 
the  actual  protection  factor  varies  from  place  to  place  inside  the 
vehicle  owing  to  the  irregular  mass  distribution  (including  shape) , 
and  to  dependence  of  the  radiation  source  intensity  on  height  above 
the  ground.  Variations  arising  from  these,  and  other  effects  as 
well,  are  at  least  10  -  15%,  experimentally,1  so  the  complexity  of 
the  actual  circumstances  has  the  effect  of  easing  the  calculational 
problem. 

To  deal  with  the  problem  of  the  extra  mass  we  have  idealized  the 
actual  vehicle  by  a  rectangular  parallelepiped  with  square  shape  and 
having  height  equal  to  that  of  the  actual  vehicle.  The  wall  thick¬ 
ness  and  the  size  of  the  square  were  fixed  by  requiring  equality 
between  the  masses  and  surface  areas  of  the  actual  vehicle  and  the 
model  vehicle.  The  equation  of  the  protection  factor  has  an 
especially  simple  approximate  form  for  this  geometry,  and  the 
results  we  have  obtained  are  in  good  agreement  with  experiment. 

Eqs.  (6.19  and  (6.17)  give  the  values  reported  in  Reference  1, 
respectively,  for  the  "survey  correlation  factor,"  which,  for  the 
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experimental  conditions  of  ground  roughness,  was  .86  times  the 
protection  factor,  and  the  calculated  results  from  our  model,  which 
assumes  a  10%  albedo  contribution  to  the  reduction  factor. 

The  results  for  the  P7,  modeled  in  this  way,  as  a  square  rectan¬ 
gular  parallelpiped,  are  given  in  Section  6  (eqs.  (6.20)and  (6.22)). 

The  outline  of  the  report  is  this.  In  Section  2  we  establish 
the  character  of  the  physical  approach  we  have  adopced,  present 
source  and  buildup  data,  motivate  our  choice  for  ground  roughness 
parameter,  and  give  the  general  expressions  of  the  protection  factor, 
the  survey  correlation  factor  and  the  RRTF.  In  Section  3  we  ex¬ 
plicate  the  geometrical  problem  of  the  reduction  factor  (unscattered 
component)  for  a  plane  slab  of  arbitrary  orientation,  introducing 
a  certain  shape  approximation  which  exploits  essential  features  of 
the  source.  In  Section  4  we  develope  the  principal  analytical 
formulas  for  the  slab  reduction  factor  and  solid  angle  factor. 

In  Section  5  we  present  the  results  of  the  M113  model  calculation 
which  neglects  the  "extra  matter"  actually  present  but  not  contained 
in  the  walls.  In  Section  6  we  present  details  of  the  parallelepiped 
model  and  give  results  for  both  the  Mil 3  and  the  P7.  In  Section  7 
we  discuss  the  results,  and  in  Section  6  we  give  a  brief  summary  of 
what  we  have  done  and  offer  a  concluding  remark  concerning 
application  of  our  method  to  computer  programing. 
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2,  PRELIMINARY  CONSIDERATIONS 
i\ .  Physical  Approach 

We  assume  the  exposure  at  a  point  inside  the  vehicle  to  be 
a  sum  of  contributions  from  all  directions,  with  chat  from  any  given 
direction  expressed  as  the  product  of  a  buildup  factor  with  the 
exponentially  attenuated  unscattered  component  deriving  from  an 
assumed  angular  dependent  source.  For  the  source  distribution  we 
assume  a  form  derived  from  curves  given  by  Spencer,2  for  air, 
above  an  infinite  plane  of  uniformly  distributed  fallout,  simulating 
effects  of  ground  roughness  by  using  an  effective  height  above  the 
ground  which  is  augmented  by  an  amount  Ah  over  the  actual  height.2 
The  normalization  of  the  angular  distribution,  t (h , x) ,  where 
x  -  cos  0,  9  the  angle  of  the  line  of  sight  from  the  detector 
position,  of  the  incident  y  ray  with  the  downward  vertical,  is  set  by 
the  usual  condition, 
r  +1 

L(h)  a|  x  -C,(h,x)  =  1,  for  h  =  3  feet  (2.1) 

To  neglect  of  albedo,  the  reduction  factor  for  detector  height 

h  ('hgff  =  h  +  Ah  if  there  is  ground  roughness)  due  to  a  slab  of 

thickness  t  and  surface  E. 

i  l 

is 

r.  =  r  (h)  =  B0 (t.)E  (h, t . )  (2.2) 

1  Li  1  Li  1 

where  B0 ( t^ )  is  the  buildup  factor  for  the  given  slab  material, 
whose  thickness  is  tj,  and  E  (t^)  is3 

E^,  (h,  t^)  =  J  dQg^, .  (n)  ■Mh,cos  8)  exp  (-p^ .  t^ ,  secB  0O )  ,  (2.3) 

where  p  is  the  inverse  attenuation  length  for  the  substance  of 
which  £  is  composed,  0O  is  the  angle  formed  by  the  incident  y-ray 
and  the  normal  to  E,  0  is  the  so-called  refraction  parameter,  which 
depends  on  the  slab  material,3  and  the  "step"  function  g£.(n)  is 
unity  if  the  line  of  sight  of  the  direction  indicated  by  Q 
intersects  E,  and  is  zero  otherwise.  In  the  same  spirit  the  re¬ 
duction  factor  inside  the  enclosure,  at  height  h  —  above  the  ground 
is  understood — is 

R(h)  =  aR  =  a?rr . (h eff) ,  (2.4) 

where  the  summation  extends  over  all  the  surfaces  of  the  enclosure. 
In  eq.  (2.4)  we  have  included  an  overall  factor  a^  to  account  for 
albedo,  i.e.  the  contribution  to  the  exposure  from  y  rays  b^ck- 
scattered  from  the  interior  faces  of  the  enclosure.  The  protection 
factor  within  the  enclosure  at  height  h  is  defined  by 
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P (h)  =  L( 3  feet) R-1  =  ft  \  (2.5) 

by  eq.  (2.1)-  it  is  the  ratio  of  the  exposure  in  the 
standard  unprotected  position  (3  feet  above  a  smooth  infinite  plane 
of  uniformly  distributed  fallout)  to  the  exposure  at  height  h 
within  the  enclosure  under  (assumed)  given  ground  roughness  con¬ 
ditions.  The  survey  correlation  factor,  S,  reported  in  Reference 
1,  is  a  similar  ratio,  with  the  measured  exposure  at  h  =  3  feet 
replacing  that  for  the  standard  unprotected  position;  so 

S(h)  =  L  (3  feet  +  Ah) ft”1  ^  L(3  feet  +  Ah) P  .  (2.6a) 

Since  L(h)  ,  defined  in  eq.  (2.1),  is  a  decreasing  function  of  h, 

S(h)  is  smaller  tl  an  P  (h)  .  We  define  the  RRTF  to  be  S  : 

RRTF  =  S"1.  (2.6b) 


B.  The  Source 

Y-rays  of 

(2.7) 

energy  E,  and 
we  have  re- 
t(h,x)  for 

h  =  6.6  feet  and  h  =  66  feet,  together  with  plots  of  eq.  (2.7)  for 
values  of  |i  such  that  the  location  of  the  peaks  of  the  former  co¬ 
incide  with  those  of  the  latter,  which  occur  at 

x  --  (i  h  =  x.,  (2.8) 

d 

and  with  normalization  of  (2.7)  chosen  so  that  the  heights  of  the 
peaks  coincide.  For  h  -  6.6  feet  the  agreement  for  x>  0  is 
extremely  good  and  for  h  =  66  feet  the  disagreement  is  of  order  10%, 
or  so.  Our  interest  is  in  height  values  appropriate  to  the  6.6  feet 
case. 

To  very  good  accuracy,  to  quite  a  bit  better  than  5%  for  heights 
of  interest  to  us,  we  are  justified  in  using  the  form  (2.7)  for 
x>0-  For  the  skyshine  part,  x<0  ,  we  assume  a  parabola  centered 
at  x  =  -1, 


l  (h,x)  *0^+02 

(x  +  l)2. 

(2.9) 

requiring  the  constants  c^  and  c 2 

l (h, -1) A (h,0)  = 

to  be  chosen  so  that 

.195 

(2.10) 

If  the  source  were  that  due  solely  to  unscattered 
energy  E,  we  would  have 

-f.  (h,x)  =  c  x  1exp(-ji  hx  ^)  ,  x>0, 

d 

t(h,x)  =  0  x<G, 

where  u  is  the  y-ray  attenuation  length  in  air  for 
c  is  a  normalization  coefficient.  In  Figures  1  and  2 
produced  two  of  the  curves  from  Reference  2,  plots  of 
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-0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.1 

FIG.  1  Exposure  angular  distribution  (x  =  cos  9)  above  cn  infinite  plane  of  radioactive  f'llout.  Solid  curve 

(taken  from  reference  2):  in  water,  scaled  to  air,  height  6.6  ft.;  1.12  hr.  fission  spectrum.  The  crosses  show 
the  function  given  in  eq.  (2.7),  with  h  =6.6  ft.,  and  a  ond  c  chosen  so  that  the  peak  of  the  correspondir.- 
curve  coincides  with  the  peak  of  the  curve  drawn  in  the  figure. 
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which  is  good  to  <5%  between  h  -  3.3  feet  and  h  =  33  feet.  We  fix 
the  relative  normalization  between  the  x  <  0  and  x>0  parts  by 
requiring 

S  =  . 110L,  (2.11) 

where  0 

S  =  d  x  t(h,x) ,  (2.12) 

J  -1 

is  the  skyshine  contribution  to  L.  Between  3  feet  and  33  feet  S 
rises  from  .093L  to  .125L,  from  curves  in  Reference  2;  interpolation 
gives  S  =  . 105L  for  h  =  6.6  feet,  while  direct  measurements  from 
Figure  1  have  given  us  eq.  (2.11). 

The  normalization  coefficient  c,  in  eq.  (2.7),  still  remains  to 
be  fixed;  it  is  a  buildup  factor,  and  it  depends  on  h.  Also, we  need 
to  choose  a  value  for  na.  First  the  latter:  the  contribution  to 
L  from  x  >  o#  assuming  the  form  in  eq.  (2.7),  is  cE^h)  where  E^(k) 
is  the  exponential  integral, 

oo 

E  (K)=y  dyy~1exp(-y)  ;  (2.13) 

we  equate  this  at  h  =  6.6  feet  to 

L(6.6  feet)  -  S(6.6  feet)  =  .890L(6.6  feet), 
taking  L(6.b  feet)  .87  from  Spencer2,  to  get 

c(6.6  feet)E^(n(6.6  feet))  =  .774.  (2.14) 

But,  with  x  =  k  locating  the  peak,  we  have  also 

c  •  (en)  1  =  >tp  (2.15) 

where  tp  is  the  peak  value  of  which  is  4.15  for  h  ■■=  6.6  feet; 
combining  these  gi ves 

kEi(k)L  .  _  t  =  .774  •  (4.15e)_1  =■-  .0686,  (2.16) 

1  6.6  feet 

v’hose  solution  is  k  =  .0207  (to  about  3%),  so  that* 

li"1  =  6.6  ft.  f  .0207  =  319  feet,  or  97m, 

li  =  .003134  ft?1  (2.17) 

3 


*  The  value  97m  seems  a  little  low  for  the  spectrum  from 
which  Spencer's  curves  are  calculated;  this  was  a  1.12 
nr.  fission  spectrum,  for  which,  using  spectrum  data 
given  in  Reference  2,  we  have  estimated  ±  123m,  to 
perhaps  15%.  a 
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If  we  use  eq.  (2.15),  taking  x  =  p.  h  and  reading  l  from 
Spencer’s  curves,  we  obtain  a  set  of  values  of  c  which^can  be  fi  t 
with  a  straight  line 

c  =  .2169  +  .835  x  ,  4  ft  £  h  5  42  ft.,  (2.18) 

while  for  h  =  3  ft.  we  get 

c (3  ft.)  -  .2172  (2.19) 

From  eqs.  (2.10),  (2.11)  and  (2.18)  we  obtain  c,  and  C2#  and  with 
this  the  complete  form  for  -t(h,x), 

!x  1exp(-pahx  *) ,  x  >  0, 

.264E1(liah)[.195  +  .805  (x+1)  ],  x<0.  (2.20) 

This  can  be  used  for  4  ft.  *h^42  ft.  even  though  (2.10)  and 
(2.11)  weaken  above  30  ft.,  for  the  whole  skyshine  contribution  is 
not  much  more  than  15  to  20%  of  the  total  below  h  =  42  ft.  When  we 
integrate  (2.20)  using  (2.19),  we  find  L(3  ft.)  =  1.000. 


The  source  function,  eq  (2.20),  is  not  quite  what  we  want, 
however.  We  need  to  shift  the  value  of  if  we  wish  to  compare 
with  experiments  of  Reference  1,  where  Co  °0  was  used  as  the  y-ray 
source.  The  Co60  spectrum  contains  two  principle  peaks,  at  1.17  MeV 
and  1.33  MeV.  if  we  idealize  this  as  a  single  peak  at  1.25  MeV,  the 
attenuation  length  in  air  is 


-1  A  -1  .  .  p 


(2.21) 


where  A/s  =  2  is  the  atomic  weight  to  number  ratio  for  air. 


-3  3 

p  =  1.312  x  10  gm/cm  is  the  density  of  air  (STP) , 

mass  of  the  proton,  a. .  the  Thompson  cross  section, 

th 


m  the 
P 


8n  2  2  2 

0  =  "y  r^,  rQ  =  e  /i^  c  =  classical  electron  radius, 

and  f(E)  is  the  Compton  cross  section  for  y  energy  E,  in  units  of4 
°th  *  We  get 


(1.25  MeV)  =  131.7  m  =  431  ft 
w  =  .002314  ft"1. 


(2.22) 


which  are  somewhat  different  from  the  values  appearing  in  (2.17). 
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To  convert  to  the  Co  spectrum,  we  simply  take  over  eq. 

(2.20),  regarding  the  right  side  as  a  function  of  n  -  h,  valid 

-1  3  -1 
between  =  .003134  ft  x  4  ft  =  .01254,  and  =  .003134  ft  x 

42  ft  =  .1316,  where  the  k  dependence  of  c  is  indicated  by  eq. 
(2.18).  There  is  an  overall  normalization  change  affecting  the 
coefficients  in  eq.  (2.18),  however,  arising  from  the  fact  that  the 
normalization  condition,  eq.  (2.1),  holds  for  any  p  .  We  have,  at 
h  =  3  ft.,  3 

c ( 3  ft,  .002314  ft'1)  =  N  •  c ( 3  ft,  .003134  ft"1)  (2.23a) 

while  from  eqs.  (2.1)  and  (2.20)  with  p  -  .002314  ft  1, 

-1  a 

c ( 3  ft,  .002314  ft  )  x  1.1222  (.006942)  -  1  (2.23b) 

whence 

N  =  .938  (2.24) 

We  asstune  that  the  same  factor  scales  the  source  to  the  new  value 
of  for  other  values  of  h  and  we  get 

c  =  .2034  +  .783u  h,  5.5  ft.  <h<  55  ft., 

cl 


=  .2038,  h  =  3  ft.  (2.25) 

The  limits  of  validity  to  (2.25)  are  determined  as  and 

*_/u  (compare  eq.  (2.18)). 

^  a 


Summarizing,  the  source  is  given  in  eqs.  (2,20)  and  (2.25) 


for  the  case  of  a  Co  source 


Buildup  Factor  for  Transmission  through  a  Wall, 


Taking  the  data  of  Reference  3,  we  can  fit  the  buildup  factor, 
as  a  function  of  (normal)  wall  thickness,  to  a  straight  line.  For 
X  out  to  2  mean  free  paths  (which  is  all  we  need)  we  have,  for 
concrete,  or  aluminum 

B  =  1  +  . 76X,  (2.26) 

o 

when  X  is  the  mass  thickness  in  units  of  p  1.  For  aluminum, 

using  eq.  (2.21)  with  p  =  2.72  gm/cm^,  we  have  |i  1  =  1  =  2.60 

in,  so  that 

B  =  1  +  .292  t  (2.27) 

o 

where  t  is  the  wall  thickness  in  inches.  This  is  good  to  a  percent 
or  so. 
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D.  Ground  Roughness 

The  data  of  Reference  1  permit  the  determination  of  the 
ground  roughness  parameter.  Ah,  which  existed  under  the  conditions 
of  the  experiment.  Figure  3  shows  a  plot  of  the  reported  unprotec¬ 
ted  exposure  as  a  function  of  height  h,  normalized  to  unity  at  3 
feet,  and  Figure  4  shows  four  theoretical  curves  of  L(h+Ah) ,  for 
Ah  =  0,3,6,  and  9  feet,  with  zeros  of  ordinate  chosen  so  that  all 
curves  pass  through  the  same  point  at  h  =  3  feet.  When  the  two 
figures  are  superposed,  the  experimental  curve  lies  on  top  of  the 
theoretical  curve  having  Ah  =  3  feet.  For  comparison  purposes 
therefore  we  adopt  this  value  for  the  ground  roughness  parameter 
throughout .  We  notice  that  if  h  <  5  feet  then  hc.002314x(5+3)  = 
,0185«1,  a  fact  we  exploit  below. 
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3.  PLANE  SLAB  HAVING  ARBITRARY  ORIENTATION 
A.  Geometrical  Preliminaries 

We  choose  cartesian  coordinates  (XQ,y  )  in  the  plane  of  the  slab, 
£,  with  the  'x  -  axis*  horizontal  and  the  'y  -  axis  forming  an 

■p 

angle  a  with  the  vertical,  where  a  =  +  I  for  the  floor  (Figure  5) . 

The  origin  cf  coordinates,  0,  is  the  intersection  of  the  perpendi¬ 
cular  to  the  x oyQ-  plane  from  tt'e  detector  position  D  =  D(0,0,z^) 

where  z^>0  is  the  distance  from  slab  to  detector;  the  detector  is 
positioned  at  height  h  above  the  horizontal  plane  of  the  ground. 

In  addition  we  employ  spherical  coordinates,  with  origin  at  D  and 
polar  axis  pointing  down.  The  co-latitude  angle  of  the  point 
P(x  y  0)  at  which  the  incident  y-ray  passes  through  Z  is  8,  by  the 

Of  Og 

definition  given  above  eq. (2.1).  Ihe  unit  vector  which  points  f rom  D to 
P  is  "n  and  the  azimuth  of  'n'  in  the  plane  parallel  to  the  ground 
(through  D)  measured  from  +  'x, counter-clockwise  from  above,  is  cp, 
the  second  polar  v triable  of  the  point**  P.  In  Figure  5  we  have 
included  also  an  auxiliary  (primed)  set  of  cartesian  axes,  having 
the  same  origin  0  and  y^  -  axis  pointing  vertically  upwards.  From 


/\  /\  / 
n  =  -  y 

o 


cos  0  +  or  cos  -  'z  sin  cp)sin  0 
o  o 


(3.1) 

(3.2a) 


/\ f  /\  •  /\ 

y  =  -  z  sin  a  +  y  cos  a, 
00  o 

/\!  /\  .  /\ 

z  =  z  cos  a  +  y  sm  a 
00  o 


(3.2b) 

(3.2c) 


we  have 

in'  =  x"  cos  cp  sin  0  -  'y'  (cos  cp  cos  6-  +  sin  0  sin  a  sin  cp) 
o  o 


+  'z'  (cos  cp  sin  a-sin  0  cos  a  sin  cp)  .  (3.3) 


♦Generally,  we  use  the  carat  to  denote  a  direction  (unit)  vector 


The  spherical  coordinate  system  is  left-handed. 
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xo»  *0 


FIG.  5  Geometry  of  slob  (X)  orientation,  detector  (D)  position,  y-roy  path  (PD)  through 
the  ground,  T  .  For  the  configuration  shown,  a>0. 
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V7e  denote  the  vector  OP  by  (T  and  pf>  by  r  so  that 
-  A  -  A  P 

rP  =  rPn  =  -  zlz0 


(3.4) 


pp =  *<A> +  yJc 
rp  =(xo +  4 +  zi)l5- 

whence,  using  eq.  (3.4) ,  we  find 

x  =  x(x  y  )  =  cos  0  =  -r  1  y  cos  a  -  r  ^z.sin  a 
o,o,  p  o  pi 

-1  .  .  -1 

tan  cp  =  -x  y  sin  a  +  x  z,  cos  a, 
o  o  o  1 


(3.5) 


(3.6 


(3.7) 

(3.8) 


frlSiSP).  =  _  _  -3 
3(xo,yo}  "  1  P  ' 

so  that  integrals  of  the  type  appearing  in  eq.  (2.3)  can  be 
performed  on  the  plane  of  the  slab, 


(3.9) 


H 


dfig^r.) 


1  f* 2  f  y2  z± 

-  *r  f  |  dxody0  •  ~ 

J*iJy  i 


(3.10) 


where  the  "left"  and  "right"  edges  of  E  are  located  at  xq  =  x^ 

x  =  x  ,  respectively,  and  the  "lower"  and  "upper’  edges  by 
O  2 

y  =  y  and  y  =  y  ,  respectively.  It  we  let 
ox  o  «  ^ 

5  =  Vl1  '  =  yoZl  '  (: 

th^n  for  (3.10)  we  may  write 

i  c  if?2  r112  2  2  -3h 

2tT J  dngz ( •  •  • =  ^7J  d?J  dri  ( ^  +  T,  +  1 )  •  •  • ,  ( - 

where  =  x.z.^  and  ti.  =  y.z.1,  i  =  12. 

l  ii  i  ii  ' 


(3.11) 


(3.12) 


We  also  need  to  express  9q  in  terms  of  the  cartesian  variables? 


from 

cos  9  =  ‘n  •  ( -/z  )  , 

o  o 

with  eq.  (3.3) , 

cos  0O  =  -  cos  0  sin  a  +  sin  9  cos  a  sin  cp. 


(3.13) 


(3.14) 
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whence  with  the  aid  of  eqs.  (3.7)  and  (3.8) 

cos  0  -  (§2  +  r\2  +  1)  (3.15) 

o 

which  also  can  be  seen  by  inspection  (Figure  3) . 

B.  Basic  Geometrical  Approximation  Scheme 

There  are  two  classes  of  slab  orientation  we  need  to  considers 

(i)  | a |  well  away  from  as  in  the  case  that  E  is  a  wall;  and 

(ii)  a  =  ±  where  E  is  a  ceiling  or  a  floor.  For  the  latter  we 
adopt  the  approximation  scheme  described  in  Reference  2,  which 
replaces  E  by  another  slab,  E0,  formed  of  suitable  combinations  of 
circles  and  annular  sectors  —  so  we  discuss  only  case  (i)  here. 

The  integrand  of  eq.  (2.3)  is  peaked  at  x  =*  k  owing  to  the  form 
of  t(h,x) .  The  intersection  of  the  plane  of  E  and  the  cone  of 
constant  x  with  vertex  at  the  detector  is  a  conic  section,  and 
t(h,x)  is  constant  along  this  curve.  We  exploit  this  fact  by  trans¬ 
forming  integrations  over  E  to  elliptical  coordinates,  u, v,  in  the 
plane  of  E,  so  chosen  that  the  curve  x  =  h,  which  (normally)  is  a 
hyperbola,  coincides  with  one  of  the  curves  of  constant  vr  v=  v^. 
From  eq.  (3.7) 

x  =  -(^2  +  r\2  +  1)  **  (ti  cos  a  +  sin  a),  (3.16) 

2  2  . 

so  the  hyperbola,  which  results  when  cos  a  >  x  is 

a'2(n  +  b) 2  -  c~2  S2  =  1  (3.17) 


where,  for  x  =  x 

2  2 

b  =  sin  a  cos  a  (cos  a  -  x  ) 

=  tan  a  +  0(x2)  ,  x  <  <  1, 

2  2.,  2,  .  2  2,-2 

a  =  x  (I-*  ) (cos  a  -  k  ) 

2  /-  2,  .  2  2,-1 

C  =  (l  —  K  ) (cos  a  -  K  ) 


(3.18) 

(3.18') 

(3.19a) 

(3.19b) 


The  angle  |  i|f  |  formed  .-»y  the  'x  -  axis  and  the  asymptotes  to 
is  given  by 

2  2  -h 

tan  )7|  =  a/c  =  k(cos  a-x  ) 

|  \lf  j  ±  x  +  0 (k3)  ,  x«l. 


(3.17) 

(3.20) 
(3.20  ' ) 


Half  the  distance  between  the  foci,  and  is 
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D  =  z^(a2  +  c2)^ 

2  lc  2  0 

-  z^(l  ~  h  )  cos  a(cosa-x^)  (3.21) 

h.  z^sec  a  +  O(h^),  k«  1  .  (3.21') 

We  introduce  the  elliptical  coordinates  as  follows.  We  let 
ri  and  r2  be  distances  of  P(x^yo)  from  P^O^bZj-D)  and  P^O.-bZj+D) 

and  define 

u=  (2D)‘1(r1+r2)  ,  v=  ( 2D)  "1  (r^)  .  (3.22) 


We  need  the  following  relations,  which  follow  from  the  above  : 

T|  =  I'Z.1  uv  -  b  (3.23) 

sec  0  =  (£2  +  t)2  +  1)  1 

o 

—  1  2  2  V 

=  Dz^  [u  -  2vqvu  +  v2  -v1(l-vo2)] 2,  (3.24) 

where 


vq  =  sin  a  and  v1  =  -h  sec  a,  (3.25) 

the  latter  being  the  value  of  v  for  which  x  =  h  holds?  also 

2  2  2  2  2  1 
x=  -  cos  a(uv  -  vov^)  [u  -  2vquv  +  v  -  v^  (1  -  v^  )]~'5,  (3.26) 


■?  »  (~  )  (u2  -  v2)[(u2  -1)  (1  -  v2)]"*5  .  (3.27) 

o (U, V)  z^ 

For  the  approximation  to  small  h,  we  drop  the  v.  terms  in  eqs. 
(3.24)  and  (3.26),  1 


sec  9q  =  sec  0Q  (u 
.  2 

x  =  -  uv  cos  a (u 
For  eq.  (3.10)  we  have 


2  % 

2v  uv  +  v  ) 
o 

2 

2v  uv  +  v 
o 


(3.241) 

(3.261) 


JdngE(n)  •  •  •  =  C2Sn  a  JJEdudv(u2  -  v2)  [  (u2  -  1)  (1  -  v2)  l”*2 

*3  / 

x  (u2  -  2v  uv  +  v2)  2  ***,  (3.28) 

o 

whare  the  cos  a  in  front  comes  from  the  scale  factor 
(D/z^) 2  •  (D/z^  (cf.eqs.  (3.24)  and  (3.27). 
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The  integrand  of  E  ,  in  eq.  (2.3),  is  proportional  to 
•t(h,x)  exp  (-a^sec  8  0q) 


where 


aE  ^E  • 


(3.29) 


We 


approximate  this  by  taking  $  ±  is 

sec  00  =  sec  0  -  (1-0)0  sec  0  tan  0  +  -  - 

n  o  n  n  n 


-aysec  0  9o  -ar  sec  0O/,  , . 

e  L  =  e  £  °(1  +  (1  -  0)0 


-agsec  0( 


sec  0  tan  0  + - )  (3.30) 

o  o  o 

(3.3o'  ) 


The  approximation  to  (3.30  )  is  reasonable  if  |a|  is  not  too 
large,  so  that  0  is  relatively  small  in  regions  where  -t(h,x) 
peaks.  For  x>0°we  have 

•Mh,x)e  SGC  9°  =  cx  ■'"expC-x  1(x  -  a^uv)  ]  ,  x>0  ,  (3.31) 

by  eqs.  (2.20),  (3. 24'),  (3.26'),  and  for  x<0  , 

2  T  ]_ 

■t(h,x)  =  [c^  +  c2(x  +  1)  ]  exp  (x  a^,  uv) ,  x<0.  (3.32) 


The  main  observation  now  is  that  x(u,v) 
dependent  in  regions  where  -t(h,x)  is  big. 
Ox(u.v)/3u)  .  2 

I  — i  =  -  cos  a  (u  -  2  sin  a 
v 

,  3  2. 

x  (v  -  sin  a  *  uv  )  , 


is  only  weakly  u  - 
To  see  this  we  evaluate 


u  v  + 


2.  -3/2 
v  ) 


(3.33) 


v 


where  the  spirit  of  the  approximation  is  the  same  as  that  used  to 
get  eqs.  (3.24/)  and  (3.26');  the  exact  expresrion  for  (3x/5u), 
vanishes  identically  at  v  =  v.  =-x  sec  a,  where  x  =  h, 

,  2 

v  =  0(h),  (dx/du)^  =  0(h  ).  For  a  =  0,  we  have 


end  for 


|  (3x/0u)v|  i  [  {~) 


+  1 


-3/2 


which,  even  for  |v|  =  .4.  satisfies  |d'x/du)v|=“  5%  since  u  j  1  by 
eq.  (3.22);  and  x(l,.4)  =  .36  ~  20h,  for  H  =  .018,  a  value  larger 
than  is  needed  for  our  purposes.  So  in  (3.31)  and  (3.32)  we  replace 
x(u,v)  with  x(l,v),  approximating  the  curves  of  constant  0,  on  E, 
by  curves  of  constant  v,  for  values  of  v  well  away  from  v^,  where 

this  result  is  exact.  The  results  are  worked  out  in  Section  4. 
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Finally  we  ;  nfr^'^e  a  certain  geometrical  approximation  Y  , 

to  the  shape  which  does  not  enjoy  complete  generality  but  which 
does  cover  a  reasonably  wide  range  of  circumstances,  those  of  the 
present  application,  in  particular.  We  simply  insert  constant, 
suitable  chosen  limits  u^,  U£,  ar ^  v^,  into  the  integral  on  the 

right  side  of  (3.28) so  that 


jang^o)  ■  •  ■  =  jV  JangSo  in) 


£f JVdv(u2  -  v2)[(u2-  1)  (1  -  v2)]^ 
i=l  “ 1  Vi 


Vl 

2  2-3/2 

x  (u  -  2  sin  a  •  uv  +  v  ) 


(3.34) 


where  the  sum  arises  from  the  fact  that  the  coordinates  u,v  only 
cover  half  the  xQy  -  plane,  e.g.  the  x  *0  half,  and  we  have  assumed 
that  x  <0<x  .  I?  0/  (x  x  )  ,  the  sum  in  (3.34)  is  replaced  by  a 

1  2  z  Cti  rv2 

difference;  i.e.,  a  single  term,  written  as  |  ^du  dv .  suffices. 

JUi  Jvx 

For  all  of  the  circumstances  of  interest  to  us  here  we  may 

choose  u ^  to  correspond  to  the  ellipse  passing  through  the  point 

(x  ,y  )  =  (x. ,0) ;  this  gives 
o  o  1 

u.  =  sec  0.  =  (x.  +  D2) ^/D  =  (1  +  §.  cos2a)^2.  (3.35) 

ill  i 


The  main  fact  making  this  choice  possible  is  chat  in  all  cases  here 
the  upper  and  lower  edges  of  at  yQ  =  y^  and  y^,  lie  between  the 

foci  of  the  coordinate  system,  at  P^  =  P  (0,  -  bz^  -  D)  and 

P  =  P  (0,  -  bz  +  D) .  In  almost  all  of  these  the  analogous  choice 
«  6  1. 


for  v  is  a  good  enough  approximation*,  that  v.  correspond  to  the 
hyperbola  passing  through  (xq,  y  )  =  (!)  y. ) .  Then,  by  the  second 

of  eqs.  (3.22),  ard  eqs.  (3.18  ')  and  (3.21'  ), 


v,  =  (y.  +  bz,)D 
ill 


-1 


=  n.  cos  a  +  sin  a. 

'i 


(3.36) 


*  The  exception  is  dealt  with  in  Section  5. 
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Figures  6  and  7  show  examples  for  two  of  the  surfaces  from 
illustrating  all  of  the  above  features.  The  choice  of  V 
by  eqs.  (3.35)  and  (3.36)  we  refer  to  as  the  centered  u 
approximation . 


the  M113, 
defined 
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Ml  13  (APC) 

Z,  (FRONT  PLATE) 


6  Centered  approximation  to  the  Ml  13  surface  Z4  .  P  and  ?2  are  the  foci  of  the  elliptical  coordinate  system  and 
the  *0  “  and  yc-axes  are  those  of  FIG.  5.  The  curves  of  constant  u  are  ellipses  and  the  curves  of  constant  v  are 
hyperbolas.  The  point  xQ  =  0,  yQ  =  y00  is  the  mid-point  of  P]  ?2‘  The  e^'Pse  u  =  °1  's  tangent  to  the  side  edges 
of  Z4,  the  hyperbola  v  =  vj  is  tangent  to  the  bottom  edge  (y0=  y]),  and  the  hyperbola  v  =  v2  is  tangent  to  the 
top  edge  (yc  =  y2). 
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4.  EXPRESSIONS  FOR  THE  E  's,  AND  ft. 


A.  Class  i.  Wall  Slabs  (|g|  Away  From  n/2) . 

It  suffices  to  treat  the  case  so  that  =  u2;  this 

is  the  symmetrical  case.  For  the  asymmetrical  case,  u^  /  u2>  the 

right  side  of  (3.34)  is  the  average  of  the  two  symmetrical  results 

for  u.  and  u_. 

1  2 

We  evaluate  (3.34)  first  with  the  "•  •  •"  replaced^by  1;  the 
resulting  expression  gives  the  solid  angle  factor  for  at  the 

detector  position,  namely  (2tt)  ^times  the  solid  angle  substended  at 

D(0,0,z  )  by  \  , 

1  L  °  Uj  v2 

.  cos  a  I*  ,  .  2  2V  r  ,  2  , ,  , ,  2, 

id  =  — - —  du  dv(u  -  v  )  [  (u  -  1)  (1  -  v  )  ] 

Lq  M  •*  1  J  Vx 


2  2-3/2 

x  (u  -  2  sin  a.uv  r  v  ) 


Using 


p  3  =  u-2v-1  y  p'  (sin  a)(V,  |v/uj<l, 
m=l  n 


where 


2  2  k 

p  =  p (sin  a)  =  (u  -  2  sin  a  •  uv  +  v  ) ^ 


(4.1) 


(4.2) 


(4.3) 


and  P  is  the  Legendre  polynomial  of  order  n,  the  prime  denoting 
differentiation  by  the  argument,  we  have 


f  du(u2  -1)  dv(l-v2)  ^[1- (v/u) 2]v 

IT  "  1  J  Vx 


2,  -Hr 


2,  -1 


•Y  P1  (sin  a) (v/u) n. 
n 


We  define 


P  “  -n.  2  P  *  n-1 Q 

=  duu  (u  -  1)  =  d@  cos  9 

n  J, 


(4.4) 


(4.5a) 


_  r  *  n-l . .  2.-H  r  *  .  -lw 

J  =  dvv  (1  -  v  )  =  |  dY  sm  Y, 

n  •)  Vl  J  •  i 


where 


v.  =  sin  Y  .  , 

i  i 


i  =  1.2, 


(4.5b) 


(4.6) 


with  ®2  given  by  eq.  (3.35),  and  interchange  orders  of  integration 
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and  summation  to  get 


•z  *  ^L-n  <si"  «)  <Vn  -  W„+3> 

o  n=l 

CD 

.  cos  a  V  ...  /. 

*  (2n  +  1  Pn  (sln  a  In+lJn+l* 

n=o 


(4.7) 


(4.8) 


in  which  we  have  used  the  identity  P^+^(x)  “  Pk_i^x^  =  (2k+l)P^(x). 

The  key  feature  of  the  expansion,  eq.  (4.9) ,  is  that  successive 
terms  diminish  rapidly  owing  to  eq.  (4.5b)  :  the  expansion  parameters 
are  the  Y^  because  the  lowest  order  term  of  the  Taylor  series  for 

J  (♦,,!_)  is  n-1(f.n-  *.n  ).  From  eqs.  (4.5),  (4.6)  and  (4.8)  we  have 
n  1  2  dL  J- 


.  cos  a  r  ~  .  3  .  .2 

_  -  r@  *  +  —  sin  a  sin  ©  .  ty 

Z  TT  i  2  2 

O 

*2 

-  ~(3  sin2a  -  l)(®2+sLn®2  cos  ®2^  *  '*'3  +  *  (4.9) 

'•’l 

The  solid  angle  factor  for  Z  can  be  computed  exactly  for 
comparison;  it  is  the  difference  of  the  solid  angle  factors  for  the 


rectangles  Z2  and  Z^,  having  common  upper  edge  at  yq  -  0  and  lower 
edges  at  (respectively)  y  =  y  and  y  =  y  (cf  Figure  6) , 

O  A  O  A. 


where,  from  eq.  (3.12) 


-  n  tan‘1^i.?2(l  +  "i2+  52» 


(4.10) 


(4.11) 


(4.12) 


For  the  case  shown  in  Figure  6,  =  55.2”,  y^  =  -72.1”, 

y  =  -20.8”,  and  x  =  49.1",  which  yields 


=  .1016 

The  parameters  of  E0  in  the  centered  approximation  are  obtained 
from  eqs.  (3.21  ')  ,  (3.35),  and  (3.36)  with  a  =  45°;  we  find  D  =  78", 
®2  =  .5620,  ^  =  -.2181,  and  =  .4556,  so  that  eq.  (4.9)  gives 
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u>  ^  .0852  +  .0203  -  .0050  *  .1005, 

Lo 

which  agrees  with  the  exact  result  for  u)^  to  less  than  2%. 

We  apply  to  E_  the  methods  used  to  obtain  eqs.  (4.8)  and  (4.9). 
With  the  help  of  eqS.  (3.26/),  (3.31)  and  (3.32)  the  analogue  of 
eq.  (4.4)  is 

E  -  C°tt  a  f  du(u2-l)  ^fdvd-v2)  ^[l-(v/u)  2]v  1 
ho  Jv 

i 

[0(-v)cx  1exp(-HX  ^)  +  0 (v) (c1+c2 (x+1) 2] 

GO 

.  exp  (-A  u)  y  p'  (sin  a)(v/u)n,  (4.13) 

h  u,  n 
n=l 


where  0(v)  denotes  the  unit  step  function,  which  vanishes  for  v<0, and 
\  s  “arvx  1  (4.14) 


In  the  approximation  discussed  near  the  end  of  Section  3, 

x(u,v)  ss  x ( 1, v)  ,  so  that  also  A^(u,v)  «  A^,(l,v)  .  In  parallel  with 

eqs.  (4.5)  we  define 


In  (v)  =  Jduu  n(u2-l)  ^exp  ( -Aj,  ( 1 ,  v)  u) 

7  rI8  n-1,  2.-^5  . 

Jn  =  J dvv  (l“v  )  li (v) > 


(4.15a) 

(4.15b) 


where  4(v)  stands  for  the  square  bracketed  expression  in  (4.13)  with 
x  =  x  ( 1 ,  v)  . 


We  write  In(v)  as  a  Taylor  series  in  v. 


00 

In(v)  = 

m=o 


-  *,>>• 


ml 


(4.16) 


substitute  this  with  eqs.  (4.15)  into  (4.13)  and  rearrange  the 
summations  to  get  the  analogue  of  eq.  (4.8), 


s,  cos  a  V  , 

-  l  J 

r=o 


f 


r+1  L 
m=o 


(2n+l)P  (sin  a). 


I 


n+1 


(r-n) 


(0) 


n 


(r-n)  1 


(4.17) 


Keeping  only  the  first  two  terms,  we  obtain  an  approximation  whose 
spirit  is  only  a  little  rougher  than  that  of  eq  (4.9), 


cos 


^[V°U)J1  +  a[X\o)  +  3  sin  a  (0))J?~  • 


(4.18) 
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The  analogous  approximation  to  uu^,  for  the  example  of  Figure  6, 
just  discussed,  agrees  with  to°a  few  percent  (1.4%).  We  note  that 

(0) 


»r(°)  =In(0) 


(4.19a) 


il1}  (0) 


n 


3 v“ 


•  'n-l'01 


=  a  tan  a  .  I  .  (0)  , 
l  n— l 

so  eq.  (4.18)  may  be  rewritten  as 
co 


(4.19b) 


EZ  * 


+  O  sin  a  "I  2 +  aY.  tan  a*I0  J2]  ^4*20^ 


If  a  =  0,  this  has  an  especially  simple  form; 
E„  »  7  1,  (0)J.,  a  =  0  . 

Xq  TT  X  X 


(4.21) 


To  complete  the  scheme  we  develope  expressions  for  J  and  J 
for  use  in  eq.  (4.20).  The  In  (0)  present  no  difficulty  because 
I  (0)  is  the  Sievert  integral;  changing  variables  in  eq.  (4.15a), 
,8a 


1.(0)  =  d6  exp  (-A  (1,0) sec  0)  =  S (A  (1,0) ,®,) , 


(4.22a) 


which  is  tabulated;5  we  note  from  eqs.  (3.26  )  and  (4.4)  that 

Aj  (1, 0)  =  a^  sec  a.  (4.22b) 


1q(0)  is  found  by  differentiating  (4.22a), 
d  s  (A,  0a )  r®2 

-  — 8  =  Jd0  sec  9  exp(-A  sec  0) 


"o 

,u2 


=  Jdu(u2-1)  ^  exp  (-Au)  , 


(4.23) 


and  this  can  be  determined  from  the  tabulation  of  S  to  adequate 
accuracy  by  numerical  differentiation.  Finally  by  examining  the 

expansions  of  In(0)  in  powers  of  for  n  =  0  and  n  =  2  we  can 

produce  the  approximation. 


I  2  (0)  «  -IQ(0)  +  2  I]/0)  +  2^0®*  eXp  ("a£  S6C  * 


(4.24) 


The  coefficients  of  the  expansions  of  both  sides  agree  up  to  and 
including  the  ®|  -  term. 
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Turning  to  J ^  and  J  ,  we  treat  the  direct  and  skyshine  parts 
separately,  writing 


j  =j<d)  +  J<! 
n  J  n  J  n 


(s) 


(4.25) 


(d) 


where  the  integration  in  j  extends  over  the  range  of  values  of  v 


for  which  x>0#  while  in 


jU 

n  ' 


x<o  always  holds.  In  the  integrand  of 


the  direct  contribution  we  extract  a  peak  factor, 

-(v  cos  a)  1  exp  [+  h(v  cos  a)  ^], 

and  expand  the  rest  in  powers  of  v.  Integrating  term  by  term  produces 
a  series  of  exponential  integrals,  E^(C^) »  having  small  arguments 
when  k«1 


where 


W  =  Ictd“ 


■*.< 


C .  =  k  sec  a  •  (-v. ) 


-1 


(4.26) 


in  which  v.  means  v  or  v  .  Using 
1  1  “ 

E2(ci}  =  ciX  exp  {_ci}  “  Ei(ci}' 


we  find,  for  v^cOcv^, 


»  c  sec  a  |E^(C1)  “  [e  ^sin  a(l-x(tan  a  +  %  cot  a))v]v  }/  (4.27a) 


and  for  v^<v2< 


(d)  rr  -C  T  3 

»  c  sec  a  \l_Ej_ (C)  ”  e  sin  a  (l-x.(tan  a  +  \  cot  a))vjv 

(4.27b) 


also,  for  v^<0<V2» 


c  sec  a  |-k  sec  a  E^(C)  +  j^e  ^  (1  -  k  tan  ajvj^  j 
and,  for  v  <  v_  <  0, 

X  ^ 


(4. 28a) 


(d) 


,va- 


J  2~'  ^  c  sec  a  {[_"K  sec  a  El^l)  +  e  ^  (I"*  tan  a^vJv 


(4.28b) 


For  the  skyshine  part  we  expand  the  integral  and  integrate:  for 


vi  <  0  <  v 


27 


NOLTR  73-137 


,  (s) 

»  (c^ 

V  .  V. 

J 

J1 

+  c2)v 

0 

and 

for  v 

'  v2  <  0’ 

,  (s) 

^S)  =o. 

J1 

=  o, 

(4.29) 


(4.30) 


In  eqs.  (4.27)  through  (4.29)  we  have  cropped  terms  of  order 
2 

h.  and,  in  comparison  with  E..(£.),  we  have  neglected  terms  of  order 

2  2  1  x 
v^  and  v2 


B.  Class  ii.  Floor  and  Ceiling  Slabs  (a  =  ±  ~) . 


The  method  we  use  in  these  cases  is  described  in  Reference  2, 
so  we  mainly  just  give  the  results  of  using  the  source  indicated 
in  eq„  (2.20) . 


Floor  (a  =  -  tt/2. 


We  suppose  the  detector  position  to  be  above  the  center  of 
a  rectangular  floor  E.  E0  is  a  shape  compounded  of  a  centered  circle 
of  radius  r,  and  two  equal,  symmetrically  placed,  annular  sectors 
having  radii  r.  and  r2  >  r^,  concentric  and  with  angle  t|l0 .  r^  is 
determined  by  the  requirement  that  the  solid  angle  factor 
subtended  at  D  by  the  circle  equal  that  of  a  concentric  square 
having  side  equal  to  the  width  of  E,  which  is  ui^'y  r2  is  determined 

the  same  way,  with  the  side  of  the  square  equal  to  the  length  of  E, 
the  solid  angle  factor  being  <jj  ,  and  ili0  is  fixed  by  requiring 


Ceiling  (a  =  +  tt/2)  . 


Since  the  skyshine  is  only  about  10%  of  the  total,  a  coarser 
approximation  will  do.  We  let  E^  be  a  circle  whose  radius,  r^,  is 

so  chosen  that  This  may  be  regarded  as  a  special  case  of 

the  surface  used“ctbove  to  represent  the  floor,  one  in  which  =  r^ , 


If  E  already  is  a  square,  u  '  =  id.  „ 
automatically  held,  and  E  becomes  a  circle. 


a’  and  r2  =  r1 
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With 

ri  *  Vi  5  Z1  tan  V  1  -  ^2,  (4.31) 

we  have 

sec  0.  .  1=1,2,  ,4.32| 

Vn"  (<V'-V)"A(^-V).  «-33) 

Expressions  for  the  uj^  may  be  produced  easily  from  eq. 

(4. 12'  .  From  the  form  (see  eqs.  (2.3),  (2.20),  (3.14),  (3.15), 

3 . 12)  )  , 

g r 

EI  =  2ild§Idr'(§2+ri2+1)~1  c  exP  [-A(?*?+n2+l),5‘1|  a=-n/  (4.34) 

5i  Th  J,  2 


(4.34) 


where 


A  =  a^+H , 


we  convert  to  integration  over  Y  ,  which 


(4.35) 


giyes 

^•2 


A.  1  \ 

r  i  2  i  _  ^2 

£  C  JdX>  (x  +1)  exP  [-W  +1)*]  +(t?o/TT)cJdXX(x2+l)"^exp(-A(X2+l)^] 


and  from  the  form, 
1  «Tla 


ES=lV|  d^dr)(§2+n2+l)"W2  {c1+c2[lr(e2+n2+l)"'sl2]. 


(4.36) 


=E1(A)-[l-(cp0/n)  ]E1(A  sec  9*  )  -  (cpj  /n)  E^  (A  sec  02  )  ;  (4.37) 


*  exp[-As(§2f112+i)!2'l  a  *  + 

derived  the  seme  way  as  eq.  (4.34),  we  find 

EE0  =[  -c1{y_1e>:p  ( +  aEE1(as.y)}-c.,{y~1exp(-aEy> 
x  [u+av+K2)  -  (1+fasiy"1+p'’2] 


y=  sec  9i 


(4.38) 


(4.39) 
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C.  Approximation  of  the  Reduction  Factor. 


The  general  approximate  form  for  R  is 


R  =- 


R 

o 


B  .  E_ 
01  £, 


Ol 


(4.40) 


In  most  cases  the  £  choices  described  in  this  section  and  the  last, 
e.g.  the  centered  approximation,  are  adequate  approximations  to  £, 
giving  results  which  are  off  only  by  a  few  percent.  But  Figure  8 
shows  an  example  in  which  the  centered  approximation  to  a  rectangular 
representation  of  a  slab  £  entails  a  sizable  distortion.  Here  we 
may  retain  eq.  (3.35),  but  change  the  values  of  the  v^  so  that  can¬ 
cellation  in  E^,  ,  of  errors  from  portions  of  £  near  the  middle  of 

and  exterior  to  £  ,  with  errors  at  the  corners,  in  £  but  exterior 
to  £  can  occur.  This  results  in  v^  slightly  smaller  in  absolute 

value  than  values  stemming  from  eq.  (3.36).  The  best  values  of  v^, 

for  example,  can  be  estimated  by  equating  the  ratio  of  the  areas  of 
the  (right)  "triangles"  with  centroids  A  and  B  in  the  figure  to  the 
values  of  the  integrand  of  E  at  these  two  points.  Tnis  condition 

gives  the  base  of  triangle  A,  which  determines  v^.  This  is  the 

method  we  have  used  to  improve  the  choice  of  £q  for  the  long  side 

walls  of  the  M113  in  the  next  section.  For  identification  purposes 
we  call  this  the  "reshaping  approximation". 
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5.  "WALLS  ONLY11  MODEL  OF  THE  M113. 

In  the  phrase  "walls  only"  we  intend  that  the  floor  and  ceiling 
surfaces  also  be  included.  Figure  8  shows  the  indexing  convention 
for  the  surfaces  of  the  M113.  To  achieve  the  idealization  implied 
by  the  contents  of  the  figure  we  have  identified  the  inner  framing 
plates  with  the  hull  walls,  ignored  a  variety  of  surface  bumps  and 
irregularities,  and  represented  the  sides  of  the  vehicle  as  flat. 
Figure  9  is  a  sketch  of  the  vehicle  cross  section  as  seen  from  the 
rear,  together  with  some  aspects  of  the  idealization  implied  by  our 
procedure. 


We  choose  the  detector  position  to  be  30  inches  above  the 
floor,  or  45.9  inches  above  the  ground,  so  that  heff  =  (36  +  45.9) 
in  =  6.8  ft,  and  *  =  .0158.  The  detector  is  centered,  equidistant 
from  the  two  sides,  and  also  more  or  less  is  centered  with  respect 
to  front  and  rear,  beneath  the  cupola  opening.  In  ‘Table  1  we  list 
the  data  for  each  of  the  2. ,  much  of  which  partially  is  determined 
by  choice  of  detector  position.  In  Table  2  we  list  the  computed 
values  of  the*  E„  ,  the  buildup  factors,  and  the  r„  . . 

’-oi  So1 


By  a  process  of  inspired  scrutiny  involving  eq.  (3.30)  we  have 
estimated  crudely  a  1.05  refraction  correction;  also,  we  assume 
a  10 ya  albedo,  so  that  from  eq.  (4.40)  and  Table  2, 


R  =  1.155R  -=  .6521 
P  =  ft-1  =  1.53 


(5.1) 

(5.2) 


*  For  Z-(Z_)  we  have  computed  E  for  the  detector  position  above 
6  7  X.q 

(below ^  the  center  of  E^(Z^)  for  simplicity.  This  is  almost 
exaxtiy  the  case  for  2  so  the  error  here  is  small;  it  is  greater 
for  Z_,  but  this  is  tolerable  because  the  total  skyshine  con¬ 
tribution  is  small,  anyway. 

Also  we  have  divided  Z^(=Z2)  into  2^  +  2*,  where  2^  is  the 

portion  of  above  the  horizontal  plane  containing  the  detector 
and  2*  the  portion  below.  For  2|  we  have  taken  u^  ={=  u2  to  take 

account  of  the  shape  irregularity.  By  the  remark  in  paragraph 

1  of  Section  4,  E  .  _  i(E  ♦  +  E„  )  where  L  and  R  refer  to  the  left 

L1  L1L  “1R 

(u  =  u^)  and  right  (u  =  u2)  portions  of 
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ib  are  big,  in  relation  to  P|?2  =  2D  =2z^  sec  a  .  The  illustration  also  shows  the  indexing 
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Ml  13  (APC) 

REAR  CROSS  SECTION 


FIG.  9  Ml  13  vehicle  cross  section,  from  the  rear,  showing  the  double  wall  structure. 

The  dotted  lines  trace  out  the  rectangular  cross  section  for  the  "walls  only"  model. 
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V1 

-.520 

0 

-.3539 

-.2164 

-.3575 

— 

— 

V2 

0 

+  „420 

+.2871 

+.4407 

-.1428 

— 

— 

®1»®2 

1.0495 

.5189 

.5620 

.4691 

— 

— 

*1 

-.547 

0 

-.3618 

-.218 

-.3656 

— 

— 

'"2 

0 

+  .433 

+.2912 

+  .456 

1433 

— 

Table  1.  Slab  data  for  computing  the  M113 

protection  factor  in  the  "walls  only" 
model.  Entries  for  dimensioned  quantities, 
in  rows  2-6,  all  are  expressed  in  inches; 
entries  in  the  last  three  rows,  10  -  12, 
are  in  radians. 
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=3 

S4 

E6 

.0516 

.0272 

.0101 

.0766 

.0144 

1.438 

1.438 

1.438 

1.438 

1.475 

.0742 

.0391 

.0145 

.1102 

.0212 

NOLTR  73-137 


From  eq.  (2.6),  using  L(6.8  ft.)  =  .86  we  find 

S  =  1.32,  (5.3) 

which  is  considerably  smaller  than  the  average  value  reported  in 
Reference  1, 

S  =  3.41,  (5.4) 

where  variations  of  detector  position  did  not  involve  deviations 
of  more  than  10  to  15%. 

The  results  of  the  simplified  model  of  Section  6  support  the 
hypothesis  that  the  discrepancy  between  (5.3)  and  (5.4)  is  due  to 
neglect  of  shielding  effects  from  material  present  additionally  to 
that  which  resides  in  the  vehicle's  surfaces. 
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6.  PARALLELI PIPED  MODEL 


the 

is 


7 

The  combat  loaded  welch t  of  the  M113  is  24,000  pounds  while 
weight  of  material  assumed  for  the  calculation  of  Section  5 


W 


shell 


kl 


(6.1) 


where  =  2.72pH  o"  2,72  X  62,4  #/ffc3  “  117  */ft3'  the  density 

of  Aluminum,  A.  is  the  ar«-  of  £.  ,  and  t.  its  thickness.  The 
i  i  i 

volume  factor  in  eq.  (6.1)  for  the  M113  is 

Z  A.t.  =  57.96  ft3  (M113)  (6.2) 

i  1  x 

so  that  W  ,  , ,  =  9838  ,  which  is  only  41%  of  the  combat  loaded 

shell 

weight  given  above.  In  light  of  this,  the  disagreement  between 
(5.3)  and  (5.4)  is  hardly  unreasonable.  We  proceed  to  a  simplified 
model  and  calculation  for  the  M113,  based  on  the  information  and 
experience  gained  so  far. 


We  idealize  the  vehicle  by  a  closed  shell  of  uniform  thickness 
and  (possibly  heterogeneous)  composition.  The  attenuation  constant 
for  the  unscattered  component  of  the  radiation  due  to  atoms  of 
species  k  is 

he  =  Vk  (6-3) 

when  is  the  number  density  for  atoms  of  species  k  and  the 
corresponding  y-ray  cross  section.  We  define  an  effective  thickness. 


‘eff  =W/«>MA)' 


(6.4) 


where  W  is  the  total  weight  of  the  fully  loaded  vehicle  and  A 
the  surface  area  of  the  shell  used  to  model  it.  Then  the  un¬ 
scattered  component  penetrating  the  wall  is  proportional  to 


is 


exp(-nte^),  where 

where 

P 

Vk’Cff 

(6.5) 

'  hM.  |  £kV 

(6.6) 

with 

fk  -  V" 

(6.7a) 

^k  (Z/AV 

(6.7b) 

37 


NOLTR  73-137 


wherein  is  the  weight  of  material  of  species  k,  Z  is  the  number 
of  electrons  on  the  indicated  atom  and  A  the  corresponding  atomic 
weight. 

The  composition  of  the  M113  was  not  available  to  us,  so  to 
estimate  p  we  have  modeled  the  material  breakdown  as  follows: 

# 

»  ^ * 1  A  AAA" 


Aluminum 

Fuel 

13/160^  dummies 
Other 


Total 


10,000 

2,000lr 

2,100# 

9,900# 

24,000# 


The 'bther" entry  includes  the  motor,  hatches,  gun,  treads,  etc.,  of 
which  we  assume  about  2/3  to  be  iron,  with  specific  gravity  =  7.86, 
and  the  rest  aluminum.  The  results,  together  with  the  values 
assumed  or  computed  for  the  f  and  the  g^,  are  qiven  in  Table  3. 

The  result  is 


S/UM  =  1’02 


so  we  just  take 

u  =  u. 


(6.8) 


as  the  attenuation  constant  for  the  shell  wall,  and 


(6.9) 


In  the  present  model  we  have  concentrated  all  the  material 
present  into  the  shell,  disregarding  effects  of  the  actual,  quite 
complicated,  distribution  of  matter.  In  the  same  spirit  we  choose 
the  simplest  shape  to  work  with,  a  parallelepiped  having  a  square 
vertical  cross  section  of  side  L,  and  a  height  H  equal  to  that  of 
the  actual  vehicle.  We  fix  L  by  equating  the  surface  area  of  the 
square  parallelepiped  to  that  of  the  best  rectangular/parallelepiped 
representation  of  the  vehicle's  actual  shape,  which  we  estimate 
from  the  Master  Diagram6: 


length  = 

162.3  in.^ 

width  = 

83.0  in.' 

>  (M113) 

(6.10) 

H 

54.3  in.J 

This  gives  A  =  372 

£t2’  ‘eft  ‘ 

4.56  in.  ,  L  =  118.1  in. 
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Species 

Sp.  Grav. 

V* 

£* 

gk 

Aluminum 

to 

• 

13,400 

.56 

1.00 

Iron 

7.86 

6,500 

.27 

.95 

Dummies 

1.8 

2,100 

.09 

1.10 

Fuel 

1.0 

2,000 

• 

o 

00 

1.23 

Table  3.  Estimated  material  breakdown  for  the  M113. 


For  all  the  walls  a  =  0,  so  we  may  use  eq.  (4.21) .  We  compute 
just  the  direct  contribution  to  R,  correcting  for  skyshine  and 

albedo  together  by  multiplying  by  1.1  (skyshine)  x  1.1  (albedo)  = 

1.21; 


ft  =  1.21  R 


(6.11) 


where 


(d) 


=  Bo(teff)(Ef 


.  y  „(dk 

+iSrWi 


(6.12) 


with  suffixes  f_  and  w  referring  to  floor  and  wall,  with  B0(teff) 
given  by  eq.  (2.27) ,  and  with 


(6.13) 


We  place  the  detector  on  the  vertical  axis  through  the  center 
of  the  square  so  all  the  E^d'are  equal;  this  also  makes  Eq  a 

circle  for  the  floor  slab.  1  For  the  wall  we  assume  the  centered 
approximation,  eqs.  (3.35)  and  (3.36).  We  assume  the  detector  to 
be  at  a  height  h^  above  the  floor.  Collecting  all  the  necessary 
formulas  specialized  to  this  case  we  have: 


p"1  =  ft.  =  1.21  R(d) 


Rld>=  Bo<W 


(E  +  4  E^d)  ) 
f  w 


Bo(teff}  =  1  +  *292  fcef f  (in) 


39 


(6.14a) 

(6.14b) 

(6.14c) 


NOLTR  73-137 


Ef  =  c(E1(A)  -  EX(A  sec  0^  ) 

(6. 14d) 

A  =  az  +  h,  ar  =  UMteff 

(6. 14e) 

H  =  .002314 (ft) _1  *  (hf  +  15.9  in  +  3  ft). 

( 6. 14f ) 

sec  0^  =  (1  -  w^)  1 

(6. 14g) 

•t  - 2n_1  tan_1 4<i +  2&'h 

(6. 14h) 

Sf  -  V(2hf) 

(6. I4i) 

E(d)=  tt_1  J(d)I  (O’# 
w  11 

( 6 . 14.j ) 

I  x (0)  =  S(aE,©2  =  45°) 

(6.14k) 

J 1  -  cE1^1),  c  =  -2034  +  .783k, 

(6.14 l) 

-  -1 
ci  =  "Kvi 

(6.14m) 

Vj^  ^  -2hf/L 

(6. 14n) 

The  sources,  respectively,  are  eqs.  (2.5)  and  (6.10),  (6.11), 

(2.27),  (4.37),  (4.35)  and  (3.29)  ,  (2.8)  and  (2.22)  (with  h  -  h  +  15.9 
in.  replaced  by  h  +  fth,  ah  =  3  ft  ground  roughness,  and  where  15.9 
is  the  height  above  ground  of  the  floor  slab),  (4.32),  (4.10)  and 

(4.12),  (3.11),  (4.21),  (3.35)  and  (4.22),  (4.27b)  and  (2.25), 

(4.26),  and  (3.36). 

We  have  not  included  a  refraction  correction  here  as  we  did  in 
Section  5(cf.  eq.  (5.1)  since  the  magnitude  of  the  effect*which  is 
the "softening"  of  the  unscattered  attenuation  at  large  90,  is  well 
overshadowed  by  that  arising  from  the  irregularities  of  the  actual 
mass  distribution.  Also,  the  neglect  of  shape  effects  implied  by 
the  assumption  of  square  vertical  cross  section  renders  such  a 
consideration  doubtful.  Finally,  a  5%  correction  can  be  absorbed 
by  modelling  the  vehicle  so  that  L  is  a  little  smaller;  for  the 
M113^  changed,  slightly  smaller,  dimensions  of  the  "best"  rectangular 
parallepiped  representation  of  the  vehicle's  actual  shape  involve 
estimates  as  reasonable  as  those  given  in  (6.10).  A  similar  remark 
holds  for  the  P7. 

We  see  from  eqs.  (6.14)  that  P  and  ft  are  fixed  as  functions  of 

h.  when  t  ...  and  L  are  given,  as,  e.g.,  in  (6.10)  for  the  M113.  The 
f  ef  f 

answers  ,  for  two  choices  of  detector  height,  hf,  above  the  floor  of 
the  vehicle,  are 
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so  that 


ft  =  ft(h  )  =  .2734 
=  .2544 


P(25.9  in)  =  3.66 
P(45.9  in.)  =  3.93 

and,  by  eq.  (2.6),  S  =  .86  P, 

S(25. 9  in. )  =  3.15 
$(45.9  in, ) =  3.38 


h^  =  25.9  in, 

=  45.9  in. 

(6.16) 

(6.17) 


! 


(6.15) 


From  Reference  1,  for  positions  18  in  from  the  (two)  side  walls, 
experiment  gives 

S(47 .25  in,)  =  3.59,  3.07  ) 

%  I  (6.18a) 

S(49. 25  in.)  =  3.76,  3.24  ) 

which  gives  an  idea  of  the  experimental  spread  arising  from  the 
matter  distribution.  The  average  of  these  values  of  S  is  3.42, 
suggesting 

S(48 . 25  inj=  3.4  ±  .3.  (6.18b) 

Experiment  also  gives,  27  in.  from  the  left  wall, 

S(31  in,)  =  3.61,  (6.18c) 

a  value  within  the  indeterminacy*  indicated  in  (6.18b). 

Comparing  these  values  with  our  computed  result,  eqs.  (6.17),  we 
see  that  the  agreement  is  good.  ("Too  good",  of  course;  a  15% 
albedo  weakens  the  agreement  by  5%.) 

g 

Turning  to  the  marine  corps  LVTP7,  we  take 

length  =  293.0  in.J 

width  =  103.6  inj  (P7)  (6.19) 

H  =  78.0  in.) 

2 

which  gives  A  =  851  ft  ,  t  ff  =  4.33  in,  L  =  181.6  in,  and 


hf  =  25.9  in.:  ft  =  .2562,  P  =  3.90,  S  =  3.35 
hf  =  45.9  in.  :  ft  =  .2471,  P  =  4.05,  S  =  3.48 
*As  distinguished  from  "uncertainty!" 


(6.20) 
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We  adopt  for  the  residual  radiation  transmission  factor  (RRTF) 
the  following: 

M113:  RRTF  =  S  =  (3.4)  =  .29  experiment 

S**1  =  (3.5)"1  =  .29  theory  (6.21) 

P7:  RRTF  =  S”1  =  (3.7) -1  =  .27  theory.  (6.22) 

Accordingly  we  recommend  the  adoption  of  the  value 

RRTF  =  .3  (LVTP7) .  (6.23) 
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7.  DISUCSSION 

The  RRTF  for  both  vehicles  is  0.3  despite  the  fact  that  the  P7 
weighs  more  than  twice  as  much  as  the  M113.  Note,  however,  that 
M113  P7 

t  ,,  =  4.56  in  and  t  =  4.33  in,  so  the  model  wall  thickness, 

eff  eff 

which  is  a  principal  determining  factor,  is  nearly  the  same  for  the 
two  cases.  This  reflects  the  fact  that  the  extra  weight  is  a  size 

effect,  LMl13  =  118.1  ia  and  hMU3 

=  54.3  in.  being  somewhat  less 

than  the  corresponding  values  for  the  P7,  Lp^  =  181.6  in.  and  Hp^  = 

78.0  in.  Moreover,  since  the  floors  of  both  vehicles  are 

15.9  in.  from  the  ground,  equal  values  of  hf  and  h  mutually  imply 

each  other;  and  since  HP7  >  H^^3  proportionately  more  of  the  total 
mass  of  the  P7  than  of  the  M113  lies  above  the  horizontal  plane 
through  the  detector  position.  This  means  that  much  of  the  P7's 
added  mass  is  positioned  so  that  it  protects  only  against  the 
relatively  ineffectual  skyshine  "threat"  (0>tt/2)  „ 
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8.  SUMMARY  AND  CONCLUDING  REMARK 

We  have  computed  the  residual  radiation  transmission  factor  for 
the  Marine  Corps  LVTP7  to  be  0.3,  based  on  a  model  of  the  mass 
distribution  whose  results  for  application  to  a  similar  vehicle,  the 
army  APC  M113,  agree  with  experiment. 

We  have  accomplished  this  by  developing  a  new  analytical 
approximation  scheme  for  transmission  through  a  rectangular 
(wall)  slab.  The  scheme  has  a  wide  range  of  applications,  to 
the  case  where  the  plane  of  the  slab  is  not  too  close  to  the 
horizontal:  specifically,  |aj  not  too  close  to  cos“lx  . 

The  method  employs  the  approximation  that  the  conics  determined  by  the 
intersection  of  the  cones  of  constant  0  with  the  plane  of  the  slab 
coincide  with  the  hyperbolic  arcs  of  a  natural  elliptical  coordinate 
system  in  the  plane  of  E.  The  coordinate  system  is  chosen  by 
requiring  the  conic  determined  by  cos  0  =  k ,  along  which  the  source 
angular  distribution  maintains  its  peak  value,  to  coincide  with  the 
corresponding  hyperbolic  arc.  Finally,  E  is  replaced  by  an  approxima¬ 
tion,  E  ,  more  convenient  for  calculation:  the  boundaries  of  Eq 
are  segments  of  ellipses  (sides)  and  hyperbolas  (top  and 
bottom)  "parallel"  to  the  curvillinear  axes  of  the  natural  coordinate 
system  for  the  problem. 

We  remark,  finally,  that  the  use  of  natural  coordinates  should 
simplify  mesh  size  problems  encountered  in  programming  the  surface 
integrals  for  numerical  integration,  at  a  considerable  savings  in 
computer  time  and  cost.  Small  mesh  sizes  can  be  chosen  for  the 
v -integration,  probing  the  rapid  variation  of  the  source  function 
through  the  peak  near  cos0  =x  ,  and  large  mesh  sizes  will  suffice 
for  the  v-integration  since  the  v-dependence  of  the  integrand  is 
relatively  weak. 
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